Watching dark solitons decay into vortex rings in a Bose-Einstein condensate 
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We have created spatial dark solitons in two-component Bose-Einstein condensates in which the 
soliton exists in one of the condensate components and the soliton nodal plane is filled with the 
second component. The filled solitons are stable for hundreds of milliseconds. The filling can be se- 
lectively removed, making the soliton more susceptible to dynamical instabilities. For a condensate 
in a spherically symmetric potential, these instabilities cause the dark soliton to decay into stable 
vortex rings. We have imaged the resulting vortex rings. 



Topological structures such as vortices and vortex rings 
have fascinated scientists and mathematicians for cen- 
turies. In a quantum fluid, vortices have quantized 
flow around a one-dimensional core where the density 
vanishes. Vortex lines, which terminate at the bound- 
aries of the quantum fluid, have been observed in su- 

~,01, and recently 
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perfluid helium and superconductors EL 
in Bose-Einstein condensates (BECs) 
vortex rings are vortices whose cores are closed loops; 
the poloidal quantized flow pattern resembles that of a 
smoke ring. Quantized vortex rings have been produced 
in superfluid helium, and were first detected in pioneering 
work by Rayfield and Reif ||. In this paper, we report 
the first direct experimental observations of vortex rings 
in BECs. 

Solitons are localized disturbances in a continuous non- 
linear medium that preserve their spatial profile due to 
a balance between the effects of dispersion and nonlin- 
earity Dark solitons have been previously created 
in single-component BECs using phase-imprinting meth- 
ods These structures are characterized by a local 
decrease in fluid density (the depth). The macroscopic 
quantum phase of the BEC differs on either side of the 
soliton; for a completely dark (black) soliton, the depth 
is 100% (a complete absence of fluid), the phase offset 
is 7r, and the soliton velocity is zero. We have created 
black solitons (having a fluid-free nodal plane) in nearly 
spherically symmetric 87 Rb BECs, and have observed the 
subsequent decay of the solitons into vortex rings. 

Unlike vortices and vortex rings, whose stability is en- 
sured by Kelvin's theorem soliton stability depends 
on the nonlinearity and geometry of the medium [|| . Op- 
tical dark solitons in self-defocusing nonlinear media have 
been observed to decay into optical vortices via a 'snake 
instability' ||, confirming predictions JTc| ] . In general, 
dark solitons in BECs are also expected to be inherently 
dynamically unstable pr|-p^| . Because both the velocity 
and depth of a soliton are determined by its phase off- 



set JTsJ , small local perturbations of either the strength of 
the nonlinearity or of the depth give the soliton a corre- 
sponding nonuniform transverse velocity profile. Once 
the soliton begins to decay, it quickly breaks up into 
more stable structures. Within our spherically symmetric 
BECs, the expected decay products are concentric vortex 
rings 0. Because the theoretical background for soli- 
ton decay in condensates has been previously described 
in detail in Ref. j^J, we limit our paper to primarily a 
phenomenological presentation of data. Evidence for soli- 
ton decay is obtained through direct images of the BEC 
density distribution. We provide a visual comparison of 
data with numerical simulations of soliton decay. 

Experimental techniques - We create dark solitons us- 
ing concepts that we have previously applied to mak- 
ing singly quantized vortices in BECs &f[(|. The 
magnetically trapped BEC can exist as a superposi- 
tion of two internal components, hyperfine states |1) = 
\F = l,m F = -1) and |2) = |F = 2,m F = 1). Conver- 
sion between the two components is achieved with a mi- 
crowave field slightly detuned from the |1) to \2) internal 
conversion energy. A "modulation beam," an off- 
resonant laser beam whose small spatial focus provides a 
requisite spatial selectivity, is dithered rapidly across the 
BEC. The laser beam induces a small ac Stark shift in the 
|1) to \2) transition energy. The phase and amplitude of 
this effective fm modulation varies from point-to-point in 
the sample in a spatial pattern determined by the time- 
varying position of the optical modulation beam; thus the 
microwave-induced inter-component conversion varies in 
phase and amplitude across the BEC. 

In our experiment, we start with a uniform-phase (i.e., 
ground state) condensate of component \2). To make a 
soliton, we choose a modulation beam pattern such that 
inter-component conversion is suppressed in the middle 
of the sample, whereas the upper and lower parts of the 
BEC undergo conversion to component |1) with an ini- 
tial phase offset of ir between the two parts. Uncon- 
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trolled variations in the timing of the modulation cause 
the angular orientation of the soliton nodal plane to 
vary randomly from one soliton to the next. For fur- 
ther details, see our earlier vortex work || and theory 
by Williams and Holland [jl6| . The soliton state contains 
about 3 x 10 5 |1) atoms, and the filling consists of about 
7 x 10 5 |2) atoms. The axial (vertical) trap frequency is 
7.6(3) Hz, the radial (horizontal) frequencies are 7.8(3) 
Hz, and the condensates have a Thomas-Fermi (TF) ra- 
dius of 28 /im. The solitons are created at temperatures 
of T = 23(6) nK, or T/T c = 0.8(1), where T c is the BEC 
phase transition temperature. 

After removing the coupling field, we observe the 
trapped excited-state condensates using nondestructive 
phase-contrast imaging. Our probe laser is tuned such 
that only component |1) is visible, and these atoms ap- 
pear bright on a dark background. As described in j|, 
a fillcd-core vortex appears as a dark hole in a bright 
atom cloud. Similarly, a soliton with a filled nodal plane 
appears as a dark band that divides the |1) BEC. This 
two-component BEC can be described as a bright soli- 
ton of component |2) trapped within the |1) dark soliton 
p7[ . The filling material of |2) atoms stabilizes the dark 
soliton against dynamical instabilities; the filled solitons 
are observed to last for at least 800 ms. 

To study the dynamical instabilities of the solitons, we 
first use resonant light pressure to selectively remove the 
|2) atoms that fill the nodal plane. The |2) atoms are 
adiabatically removed over 100 ms, and a bare dark soli- 
ton of component |1) remains. The soliton node and any 
soliton decay products are then too small to be observed 
while the condensate is held in the trap H). We there- 
fore remove the trapping potential and allow the con- 
densate to ballistically expand, causing the variations, or 
'defects,' in the condensate density distribution to also 
expand |h|. These density defects are then resolvable. 
We obtain a final near-resonance phase-contrast image of 
the expanded atom cloud, using a probe detuning of 20 
MHz. Expansion imaging has previously been used to 
detect bare vortex cores [flM . 

Numerical calculations - Before discussing experimen- 
tal data, we briefly describe results of numerical calcu- 
lations. The techniques employed have been described 
in detail in Ref. [p"4|| . The equilibrium configuration, dy- 
namics, and Bogoliubov excitation spectrum for a two- 
component condensate at low temperatures may all be 
obtained from the appropriate three-dimensional time- 
dependent Gross-Pitaevskii (GP) equation pl|] , with 
Ni = 3 x 10 5 |1) atoms in a state initially constrained to 
have odd parity along one axis, and < N2 < 7 x 10 5 
1 2) atoms in an even-parity state. The spherical trap 
frequency is chosen to be vq = 7.8 Hz, and the intra- 
and inter-species scattering lengths are (an , 022 , 012) = 
5.5(1.03,0.97,1.0) nm, respectively p2| |. The behavior 
of the BEC in full two-component simulations, in which 
we approximate a slow removal of the |2) atoms, is found 
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FIG. 1. Results of numerical simulations showing the de- 
cay of a black soliton in a BEC. The simulation corresponds 
to 3 x 10 5 87 Rb atoms in a spherically symmetric trap with 
frequency 7.8 Hz. Successive frames are shown at 50 ms in- 
tervals, with the first frame at 100 ms after the start of the 
simulation. The first row shows the density profile of the con- 
densate, integrated down an axis parallel to the soliton plane. 
The low-density regions within the cloud are also rendered 
(second row), with views perpendicular to the soliton plane. 

to be qualitatively similar to simulations where we set 
N2 = at the outset. In all cases, the soliton is found to 
undergo a snake instability, decaying into vortex rings. 

Indeed, the Bogoliubov spectrum contains modes with 
complex frequencies for all N2 considered; such modes 
have been shown to drive the soliton instability [B . For 
N2 ^ 4 x 10 5 , only one such mode remains, with imag- 
inary frequency of magnitude v ~ (4 x lO 4 /^)^- As- 
suming a soliton decay time cx v^ 1 , a filled soliton with 
N2 = 7 x 10 5 is expected to be stable for longer than 2 s, 
consistent with experimental observations. 

The results of a one-component (N2 = 0) simulation 
are shown in Fig. |} The soliton decays into three nearly 
concentric vortex rings approximately 130 ms after initial 
formation. In general, the rings tend to migrate towards 
the condensate surface, and may then grow, shrink, and 
reconnect with each other. The innermost ring is the 
most stable, remaining intact for at least 150 ms in the 
absence of a collision with another ring. Along an imag- 
ing axis parallel to the plane of a given vortex ring, the 
integrated column density reveals two dips in the den- 
sity distribution, connected by a fainter line. The ring, 
however, is not constrained to lie in any given plane, and 
may bend and tilt, producing other types of images. In 
elongated traps with aspect ratios of two or higher, sim- 
ulations indicate that a soliton can decay directly into 
nearly parallel vortex lines rather than rings; however, 
this behavior is not obtained for spherical traps. 

Experimental results - In our experimental cycle, we 
first create a filled dark soliton and obtain a non- 
destructive image of the initial trapped (filled) soliton. 
Over the next 100 ms, we remove the |2) atoms from the 
soliton nodal plane until only the |1) atoms remain. The 
trapped single-component condensate is then held for a 
variable hold time before being allowed to ballistically 
expand for 56 ms Finally, the expanded density 

distribution is recorded. 

Images of the decayed solitons typically reveal two dips 
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FIG. 2. Typical images of expanded condensates, and their 
initial states (insets) before the |2) atoms were removed. All 
images are shown at the indicated scale, (a) A vortex state, 
(b)-(e) the decay products of solitons. Image (e) was taken 
with a hold time of 500 ms; all other images were taken with 
a hold time of ms, in addition to the 100 ms |2) removal 
and the 56 ms expansion [22]. 

in the condensate density distribution, as shown in 
Fig. ||(b)-(e). While our images are consistent with the 
expected signatures of vortex rings, they are inconclusive 
evidence for the decay of solitons into vortex rings; from 
a single image containing only a pair of density dips, we 
can not differentiate between a single vortex ring and a 
pair of vortex lines. For comparison, an expanded vortex 
line is shown in Fig. ||(a). We believe that our signal- 
to-noise levels may limit our ability to see a dark line 
connecting the two clear density defects that would indi- 
cate the sides of a ring. We also note that high-frequency 
excitations visible in T — simulations would be quickly 
damped at finite temperatures, and thus not seen in the 
experiment. 

Because detection of two-dimensional topological 
structures, such as dark solitons and vortex rings, can 
be enhanced by probing along two orthogonal directions, 
we added an additional probe beam. In our improved 
apparatus, the second probe beam is orthogonal to the 
original beam, and propagates horizontally (as does the 
original). The beams intersect at the condensate. Phase- 
contrast imaging is used independently with each beam 
path, and the probe beams propagate towards different 
sides of a single charge-coupled device (CCD) camera ar- 
ray. The apparatus allows us to take simultaneous pic- 
tures of condensates from the 'front' (original) and the 
'side' (added) directions using a single camera, as illus- 
trated in Fig. ||(a). While we can observe the full depth 
of the filled soliton node along the front direction, regard- 
less of the initial orientation, the nodal plane is usually 
at an oblique angle to the side probe, and soliton con- 
trast is thus reduced. Examples of simultaneous images 
of trapped condensates from the front and side are shown 
m Figs.|(b)-(e). In the pair of images shown in Fig. ^(d), 
the filled core of a vortex is visible as the dark hole in the 
front image, and as the dark band in the side image. In 
Fig. ||(e), a filled horizontally oriented dark soliton ap- 
pears as a band across the condensate in each probing 
direction, demonstrating that the soliton splits the BEC 
into two sections. 

The use of two orthogonal probe beams confirms that 
dark solitons indeed decay into vortex rings. With the 
two beam paths, we have observed pairs of density dips in 
simultaneous side and front expanded images of decayed 
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FIG. 3. (a) An illustration of our probe beam imag- 
ing paths, (b)-(e) Simultaneous 100-^tm-square images of 
trapped condensates from the front and side imaging direc- 
tions. Shown are (b) a ground state condensate of component 
|1), (c) a shell of |1) atoms around a central ball of |2) atoms 
[21] , (d) a |1) vortex filled with |2) atoms, and (e) a hor- 
izontally oriented dark soliton of component |1), filled with 
1 2) atoms. 

solitons. In each image set shown in Figs, ^(b)-(d), pairs 
of density dips lie nearly horizontally and at correspond- 
ing locations within the two BEC images. We have also 
obtained images in which a faint line corresponding to 
a decrease in integrated column density can be observed 
between two defects in one of the expanded images, as 
shown in Figs, ^(e) and (f), providing further striking 
demonstrations of soliton decay into vortex rings. Fig- 
ures ^(e)-(g) emphasize that when the vortex ring is not 
close to horizontal, clear pairs of density dips are not seen 
in the side images. However, a weak ring-shaped halo 
may be discerned in the side image residuals of Figs. |](e) 
and (g), when the ring is in a near- vertical plane orthogo- 
nal to the side probe direction. Although the halo signal 
is only slightly above the noise of the residuals, it can be 
noticed when compared to the residuals for an expanded 
ground state condensate, shown in Fig. ^(a). Images in 
Figs. [|(b) through (f) represent the best 10% of our data 
showing solitons decaying into rings; Fig. ^(g) is typical. 

Numerical as well as experimental data show great 
variability in the types of images that can be expected 
after a soliton decays. For example, Fig. |](h) is a rare 
example showing two clear density dips in the side im- 
age, but here a clear line corresponding to the defect 
extends through the entire condensate in the front image 
and residuals. It is possible that this image shows soliton 
decay into two vortex lines rather than a ring, or breakup 
of a ring into vortex lines. We defer to a future paper 
a discussion of the diversity of our images, which may 
result from interesting vortex ring dynamics, including 
vortex ring bending, tilting, and reconnections. 

Small perturbations unavoidably induce a dark soliton 
to decay. In a spherically symmetric BEC, a soliton de- 
cays into a relatively stable vortex ring. In addition to 
demonstrating this nonlinear decay process, our exper- 
iment shows that condensates can indeed support vor- 
tex rings, which may be detected via expansion imaging. 
Future experiments may include studies of vortex ring 
stability, lifetime, and dynamics p4j, as well as investi- 
gations into vortex structures created by impurity motion 
above a superfluid critical velocity M . 
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FIG. 4. Images of expanded condensates and their initial 
states (insets), as viewed along both imaging axes. Each sub- 
figure consists of two expansion images (top row of each sub- 
figure) and two images of residuals plots (bottom row), ob- 
tained after subtracting a TF fit of each expansion image. 
The white regions in the residuals correspond to depletion of 
fluid in the BEC. The residuals are superimposed with white 
ellipses outlining the fit TF profiles of the condensates, for po- 
sition reference, (a) An expanded ground state condensate. 
(b)-(h) The decay products of solitons. The hold times (be- 
tween |2)-atom removal and BEC expansion) were ms for 
images (a)-(c),(f), and (g); 50 ms for (d) and (e); and 150 ms 
for (h). 

Possibilities to create vortex rings by other means may 
also be explored, such as by passing objects through the 
condensate [pfjf . 
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